Relaxation mechanisms of the persistent spin helix 
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We study the lifetime of the persistent spin helix in semiconductor quantum wells with equal 
Rashba- and linear Dresselhaus spin-orbit interactions. In order to address the temperature depen- 
dence of the relevant spin relaxation mechanisms we derive and solve semiclassical spin diffusion 
equations taking into account spin-dependent impurity scattering, cubic Dresselhaus spin-orbit in- 
teractions and the effect of electron-electron interactions. For the experimentally relevant regime 
we find that the lifetime of the persistent spin helix is mainly determined by the interplay of cubic 
Dresselhaus spin-orbit interaction and electron-electron interactions. We propose that even longer 
lifetimes can be achieved by generating a spatially damped spin profile instead of the persistent spin 
helix state. 

PACS numbers: 75.30.Fv, 75.40.Gb, 72.25.Dc 
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I. INTRODUCTION 

Within the field of spintronics,-~— semiconductor de- 
vices with spin-orbit coupling have attracted great at- 
tention over the past years because they offer a setting 
where electronic spin polarizations can be generated and 
manipulated in the absence of ferromagnetism or exter- 
nal magnetic fields. 4 This opens the perspective of adding 
the spin degree of freedom to the existing semiconductor 
logic in information technology without encountering the 
challenge of artificially integrating local magnetic fields in 
devices. From this applications point of view it is clearly 
desirable to maximize the spin lifetimes and coherence 
lengths in semiconductor spintronics devices. 

In this respect an ideal candidate is the persistent spin 
helix (PSH), a spin density wave state with infinite life- 
time, which exists in two-dimensional electron systems 
with Rashba and linear Dresselhaus spin-orbit interac- 
tion of equal magnitude^ due to a SU(2) symmetry of 
the corresponding Hamiltonian.— On a less abstract level 
this can be understood as the combined effect of diffu- 
sion and spin precession: the momentum-dependent spin- 
orbit field is perpendicular to the PSH wave vector, and 
its magnitude grows linearly with the projection of the 
momentum argument on the direction of this wave vector. 
If, for instance, a spin-up electron starts at the crest of z- 
spin polarization and travels at the Fermi velocity in the 
direction of the PSH wave vector, its spin precesses pre- 
cisely by an angle of 2 ir during the time it takes to cover 
the distance of one PSH wavelength. If the electron prop- 
agates off direction, the spin will still match the phase of 
the PSH everywhere because the larger traveling time to, 
e.g., the neighboring crest is exactly compensated by a 
smaller precession frequency. 

One promising progress in this direction is the recent 
realization of the persistent spin helix in a GaAs/AlGaAs 
quantum well by Koralek et ai£ They used transient 
spin grating spectroscopy^ to optically excite a sinusoidal 
profile of out-of-plane spin polarization with the "magic" 



PSH wave vector. Due to the presence of symmetry 
breaking effects in a real quantum well, instead of a state 
with infinite lifetime, two decaying modes were observed. 
Koralek et al. named these two modes the symmetry- 
reduced and -enhanced mode — the latter being the PSH. 
Although the lifetime of the observed PSH mode is not 
infinite it is still of the order of 100 ps, exceeding typical 
transient spin grating lifetimes by two orders of magni- 
tude. Intriguingly, the temperature dependence of the 
PSH lifetime displays a maximum close to 100 K. 

In order to improve the lifetimes it is necessary to figure 
out what the dominant relaxation mechanisms are. The 
temperature dependence of the PSH lifetime suggests the 
involvement of electron-electron interactions^ which are 
known to relax spin currents via the spin Coulomb drag 
effect^— However, since electron-electron interactions 
respect the SU(2) symmetry of the PSH state, they can- 
not be the sole reason for a finite lifetime but in addi- 
tion a symmetry breaking term must be present.- Here, 
we consider extrinsic spin orbit interaction 1 ^ and cubic 
Dresselhaus spin orbit interaction^ as a possible source 
of symmetry breaking as proposed by Koralek et ali^. 

It is the purpose of this work to develop a theoreti- 
cal understanding of the PSH lifetime and how the life- 
time could be enhanced. In particular we consider the 
effect of electron-electron interactions in the diffusive 
D'yakonov-Perel' regime. Regarding symmetry breaking 
mechanisms, our model (Sec. [IlJ) takes into account the 
effect of extrinsic spin-orbit coupling, which results from 
the interaction of the conduction electron spins with im- 
purities, as well as the cubic Dresselhaus spin-orbit in- 
teraction, which is known to be present in the exper- 
imental quantum well to a non-negligible amount^ In 
Sec. IIIII we derive a diffusion equation for the spin den- 
sity in our model system and discuss the contribution of 
the different symmetry breaking mechanisms. In Sec. IIVI 
we present analytical solutions for the simplified situa- 
tion where only one symmetry breaking mechanism is 
present. We propose that a spatially damped spin pro- 
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file could enhance the lifetime compared to the PSH life- 
time. For the parameters of the GaAs/AlGaAs quantum 
well used by Koralek et all (Sec. [Vj it turns out that 
electron-electron interactions in combination with cubic 
Dresselhaus spin-orbit interaction are the key ingredients 
to understand the temperature dependence of the PSH 
lifetime. Detailed conclusions and an outlook are given 
in Sec. ED 



II. MODEL 

In an envelope-function description 1 - of the conduc- 
tion band electrons in semiconductor quantum wells, the 
spin-orbit interaction takes the form of a momentum- 
dependent, in-plane effective magnetic field. The two 
dominant contributions to this field are linear in the in- 
plane momentum: The Rashba field, 1 - which has wind- 
ing number 1 in momentum space, is caused by struc- 
ture inversion asymmetry and can be tuned by changing 
the doping imbalance on both sides of the quantum well. 
The linear Dresselhaus^ contribution, in contrast, has 
winding number —1 and its physical origin is the bulk 
inversion asymmetry of the zinc-blende type quantum 
well material. It is proportional to the kinetic energy 
of the electron's out-of-plane motion and therefore de- 
creases quadratically with increasing well width. In ad- 
dition, a small cubic Dresselhaus spin-orbit interaction is 
present as well. 

Thus we write the Hamiltonian for conduction band 
electrons in the (001) grown quantum well as 



H — H + H imp + H e _ e . 



(1) 



The first term represents a two-dimensional electron gas 
(2DEG) with a quadratic dispersion and intrinsic spin- 
orbit interaction 



H ° = X ^ks> Kos's ^ks 
s,s' \k 

with the 2x2 matrix in spin space 

Ho = e k + b(k) ■ a. 



(2) 



(3) 



The ip\. s (ipks) are creation (annihilation) operators for 
electrons with momentum k and spin projection s. 
Within the standard envelope function approximation 15 
one finds = where m is the effective mass. The 
vector of Pauli matrices is denoted by cr and the in-plane 
spin-orbit field 



b(k) = b R (k) + b D {k) 



(4) 



spin-orbit interactions^ (henceforth K = 1), 
k 

bn(k) = olvf ' 



bn(k) — vf cos 20 [}' 
+ vf sin 2(f> 



lf_ /cos 36» 
7 4 ^sin36» 



sin 36* 
7 T l-cos36» 



(5) 



(6) 



Here, vf is the Fermi velocity, the angle 9 gives the di- 
rection of k with respect to the x axis and <f> denotes the 
angle between the latter and the (100) crystal axis. The 
strength of the Rashba spin-orbit field is controlled by 
a and the coefficient for linear Dresselhaus coupling /J' 
contains a momentum-dependent renormalization due to 
the presence of cubic Dresselhaus coupling, 

/3- 7 fc 2 /4, 



(7) 



where the "bare" linear Dresselhaus coefficient j3 is re- 
lated to the one for cubic Dresselhaus 7 via /3 = 7(fef ) = 
7 (ft/d) (d being the quantum well width). We assume in 
the following that the spin-orbit interaction is small com- 
pared to the Fermi energy Ep, i.e., bp/Ep -C 1, where 
&i? = b(kp) with kp being the Fermi momentum. 

Furthermore, we have included in Eq. ([T]) electron- 
impurity interactions, 



#imp = y J] ^l's' U k'ks's^ks, 



(8) 



(henceforth volume V = 1). The impurity potential is a 
matrix in spin space, 



Ukk' 



V^' {{ R l}] (l + a z ^\kxk'] 



(9) 



where the spin-dependent part arises from extrinsic spin- 
orbit interaction-^ of the conduction electrons with the 
impurity potential. In real space, the matrix oper- 
ator for electron-impurity interactions reads U(x) = 
V imp (x) + i\l/A [cr x W imp (a;)] ■ V, with V im P(a:) = 
Y^,i v i x ~ where v(x) denotes the potential of each 
single impurity, {Ri} are the impurity positions (eventu- 
ally to be averaged over) and Ao is a known material 
parameter (A = 4.7 x 10~ 10 m for GaAs). Eq. 
with V™?({Rj}) = J2jV(k' - k)e-^ k '~ k > R i, is then 
obtained by Fourier transformation. 

Finally, the Hamiltonian ((T|) contains electron-electron 
interactions, 



E 

fcl ...fc4 

81,82 



(10) 



with a Thomas- Fermi screened Coulomb potential of the 



form Vks^MM 
e(q) « 1 



■t)(|fc 3 -fci|) 



where v(q) 



Kin 



and 



contains Rashba- as well as linear and cubic Dresselhaus 



-^r with a* = ^ 47r£ 2 ° £ '- being the effective 

q a* m e z 

Bohr radius. For the GaAs dielectric constant we take a 
standard value, e r = 12.9. 
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III. SPIN DIFFUSION EQUATIONS 



electron scattering, 



A. Semiclassical kinetic equations 

Our goal is to describe the dynamics of the spin den- 
sity in real space. Using the nonequilibrium statistical 
operator method^ (see Ref. l2ll ) we derive kinetic equa- 
tions for the charge and spin components of the Wigner- 
transformed density matrix 



pk(x,t) = n k (x,t) + s k (x,t) ■ a, 



(11) 



where 



Pk: 



,(x,t) = J dre lk - r (^l,(x-r/2,t) tp s (x + r/2, *)). 



(12) 



If we restrict our calculation to the zeroth order in 
b/Ep and furthermore neglect terms that are nonlinear 
in the spin density s k (x,t)r^ the equations for charge 
and spin read 



d t n k + v ■ d x n k = J k 



imp 



2 s fc x b(k) + d t s k + v ■ d x s k = J l k v 



sre-e 
Jk J 

✓7-e-e 
" J k 



(13) 
(14) 



J£ c = 2vr J2 (2|V r i 234 | 2 - ^1234^1243) S(Ae) 
2,3,4 

[(1 - m)(l - n 2 ) n 3 m - (1 o 3, 2« 4)] , 



(17) 



J%1= 2tt£ 6(Ae) {(1-mKl-na) 



n 3 n 4 



2,3,4 
2\V 1234 \ 



n 3 1 - ni 
S3 , s A 



s 2 



-V1234V1243 , , ., 

jn 3 714 1 — n\ 1 — n 2 

(l-B-3, 2h4)}. (18) 



Here, we abbreviated j = kj (where j = 1,2,3,4 labels 
initial and final states of the two collision partners) and 
Ae = e kl + e k2 — e ka — e ki . 

In our approximation the charge kinetic equation (|13[) 
decouples from the spin kinetic equation (|14[) and is 
independently solved by the Fermi-Dirac distribution 
n k (x,t) = f(e k ) = [1 + e ^ k ' EF ^ kBT ]' 1 , where k B is 
the Boltzmann constant and T the temperature. Since 
we are not interested in charge transport or local charge 
excitations, we assume that the charge distribution is 
given by this spatially uniform solution. In the next sub- 
section we use the spin kinetic equation (fT4"]) to derive 
a drift-diffusion equation for the real space spin density, 
cf. Refs. Hllllill. 



with Vi = ki/m, where the index i — x,y labels the in- 
plane spatial directions. Note that spin and charge equa- 
tions decouple in this approximation because the gradi- 
ent terms containing d ki b(k), which would couple the 
spin and charge equations, are of higher order in b/Ep. 
Moreover, in the diffusive limit 6f t <C 1 (where r is the 
momentum relaxation time), they would yield terms of 
higher order in this small parameter bp Ti 14 ^ 24 On the 
right-hand side of Eqs. (IT3l) - (IT4l) . we have the collision 
integrals for impurity scattering, 

JT = -£ W kk , 5(Ae) An jl + ^ [(fe X k')A , 
fe' ' 

(15) 

JT = -£ W kk , 5{Ae) I As + M [k x k > ]z ( sA 

\4 (s x + s' x \ ~\ 

+^{kxkt U + s'A j, (16) 

with the transition rate W kk ' = 2 it rii \v (k' — k) | 2 , 
where n, is the impurity concentration, Ae = e k — £ k ' , 
An = rik — n k i and As = — s k i, as well as electron- 



B. Spin diffusion equations in the D'yakonov-Perel' 
regime 

In the following, we consider the D'yakonov-Perel — 
regime of strong scattering and/or weak spin-orbit inter- 
action, bp t -C 1. During the time interval t between two 
collisions which alter the momentum of an electron — and 
thereby b(k) — its spin precesses around the spin-orbit 
field only by the small angle bp t. This results in a ran- 
dom walk behavior of the spini^l In contrast to the weak 
scattering limit bpr 3> 1, the spin polarization is ac- 
tually stabilized by scattering in the strong scattering 
regime bpT -C 1: the stronger the scattering, the slower 
the D'yakonov-Perel' spin relaxation — a phenomenon of- 
ten referred to as "motional narrowing" in analogy to 
the reduction of linewidths in NMR spectroscopy due to 
disorder in the local magnetic fields. 

In the spirit of the derivation by D'yakonov and 
Perel— we will exploit the separation of the timescales 
that govern the evolution of isotropic (in momentum 
space) and anisotropic parts of the spin distribution func- 
tion. Since we deal with a spatially inhomogeneous spin 
density we also have to assume that the timescale con- 
nected to the gradient term in Eq. (Tj"4"|) is large as com- 
pared to the transport time, i.e. vp qr <C 1, where 
q is a typical wave vector of the Fourier transformed 
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spin density. Thus when speaking of "orders in bpr" 
in the following, we actually have in mind "orders in 
max{6p t, vpqT}". 

In order to solve the spin kinetic equation (|14l) we split 
off an isotropic component S(x,t) from the spin density 
Sk and expand the remaining anisotropic component in 
winding numbers and powers of momentum k, 

2 7T 

Sk = f'(e k ) S + Sfe i + Sfe i + s k .3, (19) 

m 

with 



Sk 



(20) 



n=±l 



A" 3 



5fc,i = /'(^)t3— E **»0M)e <n ', (21) 

F n=±l 

s fc ,3 = /'(<%) t£— E *fcn(as,*)e infl . (22) 

The anisotropic components of the distribution function 
arise due to the gradient term in the Boltzmann equa- 
tion and the precession around the spin-orbit field. Since 
the spin-orbit fields (0, (j^J contain terms with winding 
numbers ±1 and ±3 only these winding numbers have to 
be considered for the anisotropic part of the spin den- 
sity to lowest order in bpr. Furthermore, one can show 
that the spin density contains only the same powers of k 
as the corresponding driving terms in Hamiltonian ([3]). 
Thus we consider a k- and a /c 3 -term in the ansatz for 
the winding number ±l-terms of the spin density (|20j) 
and (|21j) . because the winding number ±l-terms of the 
kinetic equation (|14[) are the gradient term, the linear 
Rashba and Dresselhaus spin-orbit fields as well as the 



renormalization of the linear Dresselhaus term due to 
cubic Dresselhaus spin-orbit interaction. The winding 
number ±3-component of the spin density (|22|) . on the 
other hand, contains only a /c 3 -term because only the 
cubic Dresselhaus spin-orbit field contributes to winding 
number ±3 in the kinetic equation (|14|l . 

In the following we consider point-like impurities, i.e., 
isotropic scattering with t _1 = mn;v(0) 2 . Furthermore 
we assume low temperature T <C TV = Ep/ks and per- 
form a Sommerfeld expansion up to order (T/Tp) 2 in all 
momentum integrations. In this procedure we encounter 
integrals of the form (n = 2, 3, 4, 6, 8) 



de k f(e k )k n = -k F z n (T) 



(23) 



with Z2 = 1 and the Sommerfeld functions (see Ap- 
pendix [5} 



23 = 1+ YT 2+ °{n ] 

it 2 T 2 
2 T 2 

Z 6 = 1+7T —5, 



Ti 



T 2 



^ = 1 + 2^+01^ 
F \ F 



rp4 



(24) 
(25) 
(26) 
(27) 



With the goal of obtaining diffusion equations for the 
real space spin density we start by momentum integration 
of the kinetic equation, , 2 \ 2 f dk[Eq. ([13])]. using the 
ansatz (|19[) . This yields the isotropic equation for the 
isotropic component of the spin density 



d t S x = 



d t Sy 



dtS z = 



f 



2tt \2m 



kp 



7T - ( d x 8k c x + dy8k s x 



, , — f d x 5k cz + d y Sk s , 
2 7T I 2m V 



■ avp8k c>z — (5vp ^ sin 24>8k c<z + cos 2<fi8k SiZ ^J j — Z4 7 oy S x , 
avp8k s ^ z + j3vp ^sin 2<f>5k StZ — cos 20<5fc Cj2 ^ > — 24 7 ey S y , 



(28) 
(29) 



avp(8k CtX + 5k St y) + /3vp sm2(/) ^Sk CtX — Sk Siy ^ + cos 20 {^8k cy + Skg^ 



(30) 



with 



5k 
8k, 
8k 



c{s) 

(») 

c(s) 



5fc c(s ) + z 4 <5fc c(s ), (31) 
8k c ( s )-((z 4 8k c ( s ) + z 6 Sk c ( s ) - z 6 £fc c3 ( s3 )), 



where 



c 



i_3_ 

4/3 



(32) 



represents the ratio of cubic and linear Dresselhaus cou- 
pling strengths and 



Sk, 
8k 



(c3) 
(s3) 



2 Re Skuty, 
— 2 Im ^fcj^), 



2 Re Ski, 
- 21m 5fei. 



(33) 
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Eqs. (l2"8]) -(|30" |) can be seen as continuity equations 
for the spin density where the anisotropic components 
6k c r s ), <$fc c ( s )3 and 5fe c ( s ) play the role of (generalized) 
spin currents. The impurity collision integral (|16[) con- 
tains a spin-dependent part due to extrinsic spin-orbit 
interaction, which acts as a drain for in-plane spin- 
polarization with the Elliot- Yafet relaxation rate^ 

7. - (34) 

This relaxation mechanism can be understood as the net 
effect of the electron spins precessing by a small angle 
around the extrinsic spin-orbit field during the collision 
with an impurity. Since this field is perpendicular to 
the electronic motion, i.e., it points in z-direction, the z 
component of the isotropic spin density is unaffected by 
the Elliot- Yafet mechanism. 

The anisotropic components <5fe c ( s ) , Sk c r s \ and Sk c3 i s ^ 
can in turn be expressed in terms of the isotropic spin 
density Si by integrating the kinetic equation (fT4")l times 
velocity, where, this time, we omit the time derivative. 



The justification for doing so is that, in order to capture 
the slow precession-diffusion dynamics of the real space 
density, we can interpret the time derivative as a coarse- 
grained one, i.e. d t S — > AS/ At with At « bp 1 > r. 
Then the fast relaxation of the anisotropic components 
into the steady state at the beginning of each time in- 
terval At contributes only in higher order in bp t to the 
average over At. Thus, to leading order, it is sufficient to 
find the (quasi-) equilibrium solutions for the anisotropic 
coefficients. Another way of seeing this is in analogy 
with the Born-Oppenheimer approximation: similarly to 
the fast moving electrons in a molecule, which almost 
instantaneously find their equilibrium positions with re- 
spect to the slowly vibrating nuclei, the anisotropic parts 
of the spin distribution quickly adjust to the momentary 
isotropic spin density. The backaction of the anisotropic 
parts on the isotropic spin density is then well described 
using their steady state solution. 

By integrating j^Pp I dkv x{ y ) [Eq. ^SJl], equating 
terms of the same order in k and solving for the coef- 
ficients, we obtain the following anisotropic equations: 



2tt 



47r [av F (l + z 4 7swTi) - f3v F sin20(l - z 4 7swTi)] nS z -I n (d x S x + z A ^ svl Tid y S y ) 

to 



Sk C y 
5k c , z 

Sk SiV = 47r [av F (l + z 4 7swTi) + (3v F sin 20(1 - z 4 7 sw ti)] tiS z + ^£-n (d y S y + Z4% w Tid x S x ) , 
5k SiZ 

5k r . 



2 7T 

- 4:Tt/3v f ti cos 20 (1 - Z47swi"i) S z H t x (d x S y - Z47swTi<9< / S x ) , 

to 

2n 

47T (— av F + j3vp sin 20) t\S x + Aitftvp t\ cos 20 S*,, H T\d x S z , 

to 

2 7T 

- 4 7t/3w f ti cos 2(f) (1 - z 4 7 sw ti) S z -I ti (d y S x - Z4j sv ,nd x S y ) , 

TO 

2tt 



2 7T 

4tt/3vp t\ cos2(f>S x — 4ir [avp + j3vp sin 20] t\S v H 17 cL S z , 

TO 



-(5 fc„ 



^6 

4 7r/3^Csin2(/)fi(l J S w?i)S z , 



2 4 



Z6 



*c,y — Sk s ^ x — 4 7T/3w F Ccos2</)fi(l 7s W Ti)S' z 

Sk CtZ = — An /3vp(fi(sm2(j) S x + cos2cj) S y ), 
Sk SjZ = — Air (3vf^ti (cos 2(f> S x — sin 20 S^). 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 

(43) 
(44) 



The spin densities Si act as sinks and sources in the 
equations for the anisotropic coefficients <5fc±i.±3,i, 6k±i,i. 
Since the spin densities Si are determined by the initial 
conditions at t = 0, they are of zeroth order in bpr, 
whereas the anisotropic coefficients 5k±i t ±si, Sk±i^ are 
already first order in bpT. Had we included parts with 
higher winding numbers ±2, ±4, ±5, ... in our ansatz, 
these would have been generated only indirectly via the 
<5fc±i,±3,i, 5k±\ t i (all of which are already of first order in 
bp t) and would therefore be of even higher order in bp r. 



In Eqs. (I35|) - (l44l we have defined the rate of "swapping 
of the spin currents"— as 



7s- 



A k F 



(45) 



which is due to extrinsic spin-orbit interaction like the 
Elliot- Yafet rate j cy (Eq. (|34)) ), but lower order in A . 
It leads to a "swapping of spin currents" because a finite 
7 SW generates, e.g., a S y spin current in response to a gra- 
dient of the S x spin density in x direction (see Eq. (|3"9"]) ). 
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Eqs. (l35l) - (|44|) are valid to linear order in r-f sw <C 1. 

Since the anisotropic components 5k±x and 5k±x are 
related to (generalized) spin currents, the anisotropic 
equations (|3"5)) - (|4"4")) express generalized Ohm's laws. The 
effective relaxation times for the anisotropic parts of the 
spin distribution function are obtained as the inverse sum 
of the collision integrals for normal impurity scattering, 
spin-dependent impurity scattering and electron-electron 
scattering, 



1 1 

- + 7cy ze H 

7"e-e,l 
1 



7oy 



£8 



ZA r e _ e ,i 



Here, the temperature-dependent rates r r _} , , f ; 1 



(46) 
(47) 

c-c,l> 'e-e,l aC " 

count for the decay of the respective component (s& i or 
Sk,i) of the spin distribution due to two-particle Coulomb 
scattering. The rate at which winding-number-±l and 
linear-in-fc components of the spin distribution relax due 
to electron-electron interaction is 



dkxdk 2 dk 3 S(Ae) kx 



e " e ' 1 ~ ~ k B Tk F m(2ir) 4 

[1 - /(efc 3 )] t 1 - /Ofei+fea-fcs)] /0fci)/0fcJ 
{2^(1^! - fe 3 |)| 2 [cos(0 3 - 0l) fc 3 - h] 

+V(\k 1 -k 3 \)V(\k 2 -k 3 \) 
[kx + cos(0 2 - 6i) k 2 - cos(0 3 - Ox) k 3 
- cos3(0 1+2 _ 3 - 0-l) |fei + k 2 - k 3 \}} . (48) 

It is related to the spin Coulomb drag conductivity from 
Refs. lij-flll via the Drude formula. The analogous ex- 
pression for the winding-number-±l but cubic- in-k com- 
ponents reads 



dkxdk 2 dk 3 6(Ai) k\ 



f-i = I 

e ~ eA k B Tk%m(2w) 4 

[1 - [1 - f(e kl+k2 -k a )} f(e kl )f(e k2 ) 

{2\V(\k x - k 3 \)\ 2 [cos(9 3 - Ox) fcf - kf\ 

+V(\k 1 -k 3 \)V{\k 2 -k 3 \) 

[k\ + cos(6» 2 - 6»i) kl - cos(0 3 - Ox) fcf 

- cos 3(0 1+2 _ 3 - fli) |fei + ^2 - fe 3 | 3 ] } . (49) 



To find the anisotropic equations for Sk± 3 we 
follow a similar procedure as before and integrate 
J dke ±m [Eq. ([14])]. which results in 



sin 2(f) S z 

dk C 3 = 7 vp k'p 7r T3 I — cos 20 5 Z 

i cos 2cf> S y — sin 20 S x 

(cos 20 5 Z 
sm2(f)Sz 
— sin 20 — cos 20 

with 



(50) 



(51) 



T3 



1 . z s 

- + 7cy 



Z 3 23 To-o,3 



(52) 



The electron-electron scattering rate that enters the ef- 
fective relaxation time (|52|) for the winding-number-±3 
parts of the spin distribution is given by 



T c-c,3 ~~ 



k B Tk F m(2Tr)' 1 



dki_dk 2 dk 3 S(Ai) 



[1 - /(e fc3 )] [1 - /(e fel+fe2 _ fe3 )] f(e kl ) /(e fe J 

{2^(1*1 - fc 3 |)| 2 [cos3(f? 3 - 0i) k 3 3 - kf\ 

fV(|fci-fc 8 |)V(|fea-fc a |) 
[fcf + cos3(6» 2 - Ox) kl - cos 3(6*3 — 6> x ) A;| 

- cos3(0 1+2 _ 3 - Ox) \kx + k 2 - k 3 \ 3 ] } . (53) 



Finally we insert the steady-state solutions for the 
anisotropic coefficients (|35t - (|44|) and (|50|) - ([5l"j) into the 
isotropic equations (f28]) -(j30 f and obtain a closed set of 
coupled diffusion equations for the three spatial compo- 
nents of the spin density, 



<D V z - I\ 



d t S 



- 7cd Zq 

L 



■ 7ey ZA 



K xz d x - M d y 



dv z - r„ 



-K zx d x + M z d v 



-K 8 

J v zy u y 



7cd Z 6 - 7cy z A 

M z d x 



K 8 



S. (54) 



DV - T x - T y - 2 7 cd z 6 - T s 



On its diagonal the matrix operator contains the pure relaxation rate 7 ey due to extrinsic spin-orbit interaction, 
diffusion terms with V 2 = d x + dy and the Elliot- Yafet In addition, it contains the D'yakonov'-Perel' relaxation 
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rates ^x(y) and 7cd which reflect the randomization of the 
spin orientation due to precession (between the collisions) 
around the winding-number-±f and winding-number-±3 
spin-orbit fields, respectively. The S x component is re- 
laxed as a consequence of precession about the y com- 
ponent of the spin-orbit field only, and vice versa. In 
contrast, the S z component is relaxed by the precession 
about the full spin-orbit field. Thus the relaxation rate of 
S z due to precession is the sum of the ones for S x and S y , 



plus a correction T sw for processes that involves the swap- 
ping of the spin currents due to extrinsic spin-orbit inter- 
action. Due to precession there are also off-diagonal rates 
L, which couple the in-plane spin components, as well 
as several off-diagonal mixed diffusion-precession rates, 
which are accompanied by partial derivatives. 

In terms of the parameters of our model and previously 
defined quantities, the coefficients in the spin diffusion 
equation (154)) are given by: 



7cd 



K 



K 



xz(yz)((P 
zx(zy){4> 

M((f> 
MM 

m 



v 2 F 7 2 fc|r 3 , 



4 

1 9 

1 sw — — % 7sw 



8 2 



sin 2(f) H - 

or 



D-(z 4 (D + D)+( 2 z 6 D 



2 



DT l z A -'—[DT 1 z A -CDf l z & -CDTizi + C,DTi — 



24 



qo \ D T 
qo [Dt 



cos 2(f) q 
cos 2(f) qo 



D- l -Q Zi (D + D) 
D-\tz 4 (D + D) 
D-^z 4 (D + D) 



sin 2(f) 



1 



7sw qo Ti D z 4 ± 



Ti D z 4 - C h D Zq 



1 



sin 2(f)) + - -f sw n q D z 4 



1 ± - (1 - C^4)sin20 
a 



- o 7sw qo cos 2(f) - 
2 a 



Ti D Z4 - ( h D z 6 



D--(z 4 (D + D) 



1 p , 

- n 7sw n qo D Zi COS 2(f) — (I - C z 4 ) , 

2 a 



COS 2(f) — On — 

2 y ° a 



D--Cz i (D + D) 



(55) 
(56) 

(57) 

>urln). (58) 
(59) 
(60) 

(61) 
(62) 



with the PSH wave vector 

qo = Avpma (63) 
and the effective diffusion constants 

D = ~4n, b = ivlh- (64) 

At T — 0, we have z n — 1 and electron-electron inter- 
actions are absent, such that D = D. Then, if we leave 
out extrinsic spin-orbit interaction in the spin diffusion 
equation , it agrees with the one presented in Ref. [ID 
(except for the sign of L). If we further omit cubic Dres- 
selhaus spin-orbit interaction in our diffusion equation, 
it also concurs with the one of Ref. |6| provided that the 
spin-charge coupling is negligible. 

IV. PERSISTENT SPIN HELIX IN THE 
PRESENCE OF SYMMETRY BREAKING 
MECHANISMS 

In this section, we use the spin diffusion equation ([54]) 
to calculate the lifetime of the persistent spin helix in the 



presence of symmetry breaking mechanisms. We consider 
extrinsic spin-orbit interaction, cubic Dresselhaus spin- 
orbit interaction or simple spin-flip scattering as possible 
symmetry breaking mechanisms. In order to allow for 
simple analytical solutions we discuss each of the symme- 
try breaking mechanisms separately. In the case of cubic 
Dresselhaus spin-orbit interaction we neglect at first the 
renormalization of the linear Dresselhaus spin-orbit in- 
teraction (see Eq. (JT))). This is formally achieved by set- 
ting ( = in Eqs. (|5"6"1) - (|6"2"|) while keeping the 7 c d term in 
Eq. (|5"4"|) . However, we will include the renormalization of 
the linear Dresselhaus spin-orbit interaction when we dis- 
cuss a possible stationary solution and when we compare 
to the experimental results in a GaAs/AlGaAs quantum 
well in Sec. 



We choose our coordinate system such that the x axis 
points into the (llO)-crystal direction, corresponding to 
(f> = j in Eqs. (l56"l) -([6"2" j) . Considering an initial spin 
polarization, which is uniform in x-direction, then due to 
L(j) = M(j) =0 the S x component decouples from the 
S y and S z components and we can set S x = 0. For a = /3 
Eq. (1541) reduces for the remaining S y and S z components 
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TABLE I: Specification of X and N in Eq. ([65 





simple spin flips 


extr. spin-orbit int. 


cubic Dress. 


X 








N 


1 





2 



to 



d t S 



Dd 2 v -q 2 D-X 
-2 g D d y 



2q Dd y . 



(65) 



where the relaxation rates due to the respective 
symmetry-breaking mechanism are represented by X and 
an integer N according to Table Q] 

For the SU(2) symmetric situation X = there exists 
a steady state solution with wave vector q$. This is the 
persistent spin helix state. More precisely for an initial 
spin polarization of the form 



S(x,t = 0) = S* o (0, 0, cosg y) 



(66) 



similar to the experimental set-upj£ one finds that the 
time-dependent solution to Eq. (|65)) is 



So 



1] sin q y 



s x -°(y,t) = ^ 2 , (67) 

2 \[e 4 «o ut + 1] coaqoyj 

For t — > oo, i.e., in the stationary limit, this reduces to 
the persistent spin helix state. 

In the presence of symmetry breaking mechanisms, i.e. 
for X ^ 0, one can still find a steady state solution of 
the form 



S y (y) = -ye-^d singxy, 



S'o 



e v/llL (C 2 sin q x y + cos q x y) 



(68) 
(69) 



This solution is a spatially damped persistent spin helix 
state with coefficients given by 



Ix 1 



qx 



2 
(Jo 



C 2 = 



|V2S+(7V + l)e-2, 



V2S-(7V + l)e + 2, 

4y/2S- (I + AQe + 2 
e (-8(N 2 — 1) + (N — 1)3^) 

4+(3N+l)£-N(N - 1)£ 2 
-(4+(JV-1)0h" 

8 — (N — 1) 2 £ 2 - 4 (2 5 - (JV + 1) Q 
(iV-f)eV8(^+l)^-(^-l) 2 e 2: 



(70) 
(71) 



(72) 
(73) 



where f = AT/(g 2 D) and S = ^(1 + 0(1 + ^0- In the 
absence of symmetry breaking mechanisms (£ — > 0) the 



t — > oo asymptotics of Eq. (I67)) . i.e., the truely persis- 
tent spin helix state, is recovered. The spatially damped 
persistent spin helix state (|68 ]) - (|69l) could in principle be 
excited with the initial spin polarization profile 



S(x,t = 0) = S e- y / lx (0, 0, cosq x y) 



(74) 



Although the spatially damped persistent spin helix is 
clearly a steady state solution when the symmetry break- 
ing is caused by simple spin flips or extrinsic spin orbit 
interaction, it is not obvious that this applies also to the 
case of cubic Dresselhaus spin orbit interaction, since we 
have neglected the renormalization of the linear Dressel- 
haus spin orbit interaction (C 7^ 0), which might lead to 
a finite lifetime of the spatially damped state. Neverthe- 
less, even when the renormalization of the linear Dres- 
selhaus spin orbit interaction is taken into account one 
can still find a steady state solution of the form (| 6 8 1) - ([73"]) 
when the ratio of the linear Rashba and Dresselhaus spin 
orbit interactions is given by 



D 



D-Uza (D + D) 



(75) 



Then the spin diffusion equation ([Ml) can still be cast 
into the form of Eq. (|65l) when the symmetry break- 
ing rate X is redefined as X = X + q^D F(T) with the 
temperature dependent dimensionless function F(T) = 

— 1 ) . For this symmetry 



'-(- 

4 \D- 



D 2 -(z 4 D (D+D)+Cz 6 DD 



4\D2-(z 4 D(D+D) + ±Czi{D+D) 2 

breaking rate X and spin orbit couplings of ([75]) the spa- 
tially damped spin profile of the form ([68l) - ([73)l is again 
infinitely long-lived. 

This stationary state should in principle be realizable 
in the GaAs/AlGaAs quantum well used in Ref. [jj be- 
cause there the ratio of /3/a almost fulfills relation (f75|) 
at a temperature of T = 100 K. For theparameters of the 
GaAs/AlGaAs quantum well of Ref. @ the steady state 
solution ([68 |) -([73 |) would be characterized by a wavevec- 
tor of gjjr f=a go and a damping length of a bit more than 
a PSH wavelength, k, 1.06^. Although a spin grat- 
ing with such a strong spatial damping might be difficult 
to realize, it should be noted that the required damp- 
ing length is oc £ _1 j so that one can expect much longer 
damping lengths for thinner quantum well. 

We now want to consider the conventional PSH solu- 
tion. When we stick to an initial spin polarization with 
the form of a plane wave (l66l) similar to the experimental 
set-up^ the time dependent solution is given by a double 
exponential decay, 



S y {y,t) 
S z (y,t) 



So . 

"7T sm q y 



So 

y cos g ?/ 



4 gp D T « 



y/{4q 2 D) 2 + {N~l) 2 X< 



(76) 



(N - l)X[e 



v/(4g 2 #) 2 + (Ar-l) 2 A^ 



(77) 
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with the symmetry-en/ianced and -reduced lifetimes 



'E(R) 



2q 2 D + -(N 



DX 



t~V(4% 2 £>) 2 + (^-1) 2 * 5 



(78) 



In the absence of any symmetry-breaking relaxation 
mechanism, i.e., for X = 0, the proper persistent spin 
helix state is recovered (te = oo). Expanding Eq. (|78|) 
for small X/ (4 Oq D) < 1 we obtain 



TR 



(N + l) 



x- 



(N 



N -I 



N- 
1)X 



1 



AqlD 2{AqlD) 



(79) 
(80) 



The reduced lifetime tr is not very sensitive to details of 
the symmetry-breaking mechanism as long as it is weak. 
Correspondingly, the temperature dependence of the re- 
duced lifetime tr is almost independent of the symme- 
try breaking mechanism (and is given by the electron- 
electron relaxation rate t~1 1 contained in D via n , 
see Eq. (|46p V The temperature dependence of the en- 
hanced lifetime Te, by contrast, depends crucially on the 
symmetry breaking mechanism under consideration and 
thus offers a way to discriminate between the different 
symmetry breaking mechanisms. For small symmetry 
breaking terms the enhanced lifetime te is proportional 
to the respective scattering rate X^ 1 . Therefore also the 
temperature dependence of te is determined by the re- 
spective scattering rate. For simple spin-flip scattering 
X = t^ 1 we expect a temperature independent lifetime 
te due to constant r s f. For extrinsic spin-orbit interac- 
tions, X = 7 cy , to leading order in A/(4oqZ)) the only 
temperature dependence comes from the Sommerfeld cor- 
rections. Thus te decreases quadratically with temper- 
ature. For cubic Dresselhaus spin-orbit interaction one 
finds 



TE 



■ %d Z 6 1 



and therefore te is proportional to r 3 1 (see Eq. ([55 



Since T3 decreases with temperature because of en- 
hanced electron-electron scattering r~g 3 (see Eq. (|52jl) 
the lifetime te increases initially with temperature due 
to the motional narrowing effect in the D'yakonov-Perel' 
regime. The presence of the Sommerfeld function zq, on 
the other hand, leads to a decrease of te with increas- 
ing temperature. Thus for cubic Dresselhaus spin-orbit 
interaction we find that the temperature dependence is 
governed by a competition between increasing and de- 
creasing contributions. We will compare this theoretical 
interpretation with experimental results for the persis- 
tent spin helix in GaAs/AlGaAs quantum wellai in the 
next section. 
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FIG. 1: (a) Temperature-dependent relaxation rates due to 
electron-electron interactions (solid line), f~* 1 (dot- 

dashed line) and T~^ 3 (dashed line), as numerically computed 
using the experimental parameters of Ref[7|. In order to con- 
tinuously interpolate between the data points, we made a fit 
to the functional form A T 2 +B T 2 In T, which has been shown 
to be correct for the spin Coulomb drag conductivity at low 
temperatures in Ref. [ill . For comparison we also show (blue 
dotted line) the inverse transport time r~ 1 (100 K). (b) The 
resulting effective relaxation rates Tj - (solid), Tg" 1 (dashed) 
and ff 1 (dot-dashed), cf. Eqs. (g6j)-(|4TJ) and (52]). 



V. PERSISTENT SPIN HELIX IN GaAs/AlGaAs 
QUANTUM WELLS 

In order to address the lifetime of the PSH observed 
experimentally in GaAs/AlGaAs quantum wells^ we con- 
sider cubic Dresselhaus alongside with extrinsic spin- 
orbit interaction as possible symmetry breaking mech- 
anisms. We also include also the renormalization of the 
linear Dresselhaus coupling constant due to cubic Dres- 
selhaus spin-orbit interaction (( ^ in Eqs. ([5u]) -(|B2"] ) ). 
Analogously to the previous section we can set S x = 
and then the spin diffusion equation (|54l) reduces for the 
remaining components S y and S z to: 



d t S 



Y 

-K zy (7T/4)dy 



K. 



yz {ir/4)d y 

z 



Ddl 



with 



Z = T :c (7r/4) + r y (^/4) + 2 7 cd z 6 + r s 



(82) 



(83) 
(84) 



For an initial spin polarization of the form S(x, t = 0) = 
So (0, 0, cosqoy) the time dependent part of the solution 
is given by a double exponential decay, i.e., a sum of two 
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exponentially decaying terms with a symmetry enhanced 
relaxation rate te and a symmetry reduced relaxation 
rate tr given by 



(85) 



4q$K yz (ir/4)K zy (w/4). 



In order to compare our theory with the experiment of 
Ref. we need to calculate the coefficients that occur in 
Eq. (|85l) — in particular the temperature-dependent rates 



for electron-electron scattering. Fig. [TJ (a) shows Tg.^, 



fr 1 



e _ e 1 and T e _ e 3 , evaluated from Eqs. 



53]) by Monte 

Carlo integration for parameters of Ref. 7. With these 
electron-electron scattering rates we find for the effective 
scattering rates in Eqs. (|46)) - (|47|) and (|52|) the results 
depicted in Fig. QTb). 

In Fig. [3J we show the numerical results for the tem- 
perature dependence of the symmetry-enhanced and re- 
duced lifetimes te(r) where we use the experimental pa- 
rameters of Ref. 0- In particular, we take a = 0.0013 for 
the Rashba spin-orbit interaction and jvf — 5.0 eV A 
for the cubic Dresselhaus spin-orbit interaction. We 
adjust the linear Dresselhaus spin-orbit interaction to 
/3 = 1.29 a in order to maximize te for T — 75 K — the 
temperature at which also in the experiment the spin- 
orbit interaction was tuned in order to maximize te- 

At intermediate temperatures around 100 K, i.e., in 
the temperature range where our theory should be most 
applicable, we find very good agreement between our the- 
ory (solid lines) and the experimental lifetimes (dots) , see 
Fig. [2] (a). We observe a maximum in te roughly where 
the experimental points exhibit one. Also the size of te 
as well as of tr is very close to the experimental val- 
ues. Since the scattering rates due to extrinsic spin-orbit 
interaction are very small in the GaAs/AlGaAs quan- 
tum well under consideration, i.e., 7 cy /7 c d ~ 10~ 4 and 
■7"7sw ~ 3 x 10 -3 , effects of extrinsic spin-orbit interac- 
tion turn out to be negligible, see Fig. [U (b). A calcula- 
tion which includes extrinsic spin-orbit interactions and 
electron-electron interactions but excludes cubic Dressel- 
haus spin-orbit interaction (not depicted in Fig. 2) would 
yield enhanced lifetimes that exceed the experimental 
ones by a factor ~ 10 3 . 

Interestingly, the simple result (151]) for the enhanced 
lifetime, where we neglected the renormalization of the 
linear Dresselhaus spin-orbit interaction due to cubic 
Dresselhaus spin-orbit interaction, is a fairly good ap- 
proximation (see dashed curve in Fig. [5] (a)). Thus the 
simple interpretation of the temperature dependence of 
te can also be extended to the present situation. The 
formation of the maximum in te at intermediate tem- 
peratures around 100 K is caused by the competition be- 
tween two effects: on the one hand te increases with tem- 
perature due to increasing electron-electron scattering, 
which leads in the presence of symmetry breaking cubic 
Dresselhaus spin-orbit interaction to the usual motional- 
narrowing effect in the D'yakonov-Perel' regime. On the 
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FIG. 2: (Color online) (a) Temperature-dependent lifetimes 
of the enhanced (red/grey) and reduced (blue/black) modes. 
The points are experimental data from Ref. 0. Solid lines are 
the respective theoretical curves including extrinsic and cubic 
Dresselhaus spin-orbit interactions as well as electron-electron 
interactions; the thin dashed line is the simplified result from 
Eq. (|8ip. In panel (b) we zoom in on the theoretical curve 
of panel (a) close to the maximum of te- The dashed line 
is the theoretical curve without extrinsic spin-orbit interac- 
tion. Panel (c) depicts the results of a calculation, where we 
include extrinsic and cubic Dresselhaus spin-orbit interaction 
but exclude electron-electron interactions. (Also here, a com- 
parison as in (b) would show that the influence of the extrinsic 
spin-orbit interactions is marginal.) 



other hand the magnitude of Sommerfeld corrections in- 
creases with temperature reducing the lifetime te in two 
ways: (i) by increasing the effective cubic Dresselhaus 
scattering rate 7 c d z$ and (ii) by increasing the linear 
renormalization of the Dresselhaus spin-orbit interaction, 
which leads to a detuning of the Rashba and the effective 
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linear Dresselhaus spin-orbit interactions. 

The important effect of electron-electron interaction 
for the temperature dependence of the lifetimes te and 
tr can also be deduced from Fig.[2][c), where we show the 
lifetimes excluding the effect of electron-electron interac- 
tions. Obviously the initial increase of the lifetimes with 
temperature is absent for both te and tr in the absence 
of electron-electron interaction. 

At low temperatures and at high temperatures devia- 
tions between our theory and the experimental lifetimes 
are more pronounced. We suppose that at high tem- 
peratures symmetry breaking mechanisms that are not 
captured by our model (e.g. effects involving phonons) 
could become important. Furthermore, since the Fermi 
temperature in the GaAs/AlGaAs quantum well under 
consideration is only Tp = 400 K we cannot expect our 
calculation, which is based on a low-order Sommerfeld 
expansion, to be as accurate in the high temperature 
range above 200 K. The disagreement at low tempera- 
tures, on the other hand, results most likely from the 
fact that we do not take into account the temperature 
dependence of the transport lifetime but rather use the 
experimental 100 K-transport lifetime r(100 K) = 1 ps 
at all temperatures. In reality, however, the transport 
lifetime increases strongly with decreasing temperature 7 
such that bp t± > 1 for low temperatures, i.e., the system 
is outside the D'yakonov-Perel' regime and our theory is 
no longer applicable. In this low temperature regime 
other approaches which account for strong spin-orbit in- 
teraction could be usedi 29 ' 30 



VI. CONCLUSIONS 

Using a spin-coherent Boltzmann-type approach we 
have derived semiclassical spin-diffusion equations for a 
two-dimensional electron gas with Rashba and Dressel- 
haus spin-orbit interactions including the effect of cubic 
Dresselhaus and extrinsic spin orbit interactions as well 
as the influence of electron-electron interactions. Based 
on this approach we have analyzed the role of electron- 
electron interaction in generating a finite lifetime of the 
persistent spin helix state. 

Our calculation shows that the Hamiltonian has to con- 
tain SU(2)-breaking terms such as cubic Dresselhaus or 
extrinsic spin orbit interactions in addition to electron- 
electron interactions. Otherwise the persistent spin helix 
remains infinitely long-lived. We find that in this respect 
the effect of extrinsic spin-orbit interaction is negligible 
in the quantum wells used in the experiment by Koralek 
et alZ Instead, the experimentally observed temperature 
dependence of the lifetime of the persistent spin helix, 
which displays a maximum at intermediate temperatures 
close to 100 K, is caused by the interplay of cubic Dressel- 
haus spin-orbit interaction and electron-electron interac- 
tions. The formation of the maximum can be understood 
as follows: due to electron-electron interactions the scat- 
tering rates of the winding number ±3 components of the 



spin distribution function grow with increasing tempera- 
ture. Since the inverse of these rates enters the effective 
scattering rate in the D'yakonov-Perel' regime, electron- 
electron interactions increase the PSH lifetime with tem- 
perature. On the other hand, Sommerfeld corrections of 
the cubic Dresselhaus spin orbit interaction enter directly 
into the expressions for the effective scattering rates and 
thus decrease the lifetime of the PSH state with increas- 
ing temperature. Also temperature-induced deviations 
from the SU(2) point due to a renormalization of the lin- 
ear Dresselhaus coupling constant by cubic Dresselhaus 
spin-orbit interaction increase with temperature and thus 
effectively reduce the lifetime of the PSH state. Since 
these corrections due to cubic Dresselhaus spin-orbit in- 
teraction dominate for larger temperatures, whereas the 
effect of electron-electron interaction prevails for lower 
temperatures, a maximum of the PSH lifetime emerges 
at intermediate temperatures. 

Our theory reproduces qualitatively the lifetime of the 
PSH state in the whole temperature range accessed ex- 
perimentally by Koralek et alZ. For intermediate tem- 
peratures close to the maximum, i.e., in the regime where 
our diffusive theory should be valid, we find also quanti- 
tative agreement with the experimental data. 

In order to maximize the lifetime, we propose to use 
a spatially damped sinusoidal spin profile as an initial 
condition for a transient spin grating spectroscopy ex- 
periment. When cubic Dresselhaus spin orbit interac- 
tion represents the only SU(2) symmetry breaking ele- 
ment, the proposed spin density profile is infinitely long 
lived similar to the PSH state in the absence of symmetry 
breaking terms. 

It may be interesting to include also disorder in 
the local Rashba spin-orbit coupling or spin-dependent 
electron-electron scattering in order to apply our theory 
to situations where the cubic Dresselhaus spin-orbit in- 
teraction is less dominant. These relaxation mechanisms 
are currently discussed in the context of spin relaxation 
in (110) grown GaAs quantum wells i 31 ' 32 
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Appendix A: Sommerfeld functions 

From the standard Sommerfeld technique in the theory 
of the Fermi gas it is well known that the approximation 

J deg{e)f(e) = J de g(e) + — (k B T) 2 g> (E F ) 

+ 0(T 4 /T F ) (Al) 

holds, where /(e) is the Fermi distribution and g(e) is a 
function of the energy that varies slowly for e ~ Ep. In 
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the derivation of the spin diffusion equations we have to 
deal with powers of momentum k 2 ,k 3 ,k 4 ,k 6 ,k 8 . Since 
the dispersion is quadratic and the 2d DOS is constant, 
the problem reduces to (n = 1, |, 2, 3, 4) 



/>oo />oo 

/ dee n f'(e) = - dene*- 1 /(e) 
Jo Jo 





-{Ef, 
+■ 0(T 4 /T 4 ) 



. 7T 2 /fc S T 



(A2) 



Thus, fc 2 terms do not obtain any T-dependent correc- 
tions, whereas the higher powers, k 3 ,k 4 ,k 6 and k s , are 

- kp but acquire correc- 
z 3 , . . . z 8 , see Eqs. I[2"4"|)- 



-P ■ 

H>p, 



not simply replaced by 
tions in the form of the factors 
(123. 
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